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ABSTRACT 


AfflDonls has been detected In the clrcumstel lar envelope of IRC 410216 
by means of 3 Infrared absorption lines In the Vj band around 950 cm 
The lines are fully resolved at a resolution of 0.22 km/sec and Indicate 
that most of the circumstellar gas Is accelerated to expansion velocities 
around 14 km/sec within a few stellar radii. The NH^ profiles Indicate 
a rotational temperature between 400 and 700 K, an density between 
lO^cm ^ and lO^^cm and an NH^ column density of lO^^cm”^. The 
density Indicates that the mass of the circumstellar envelope within a 
1 arcspc '•adlus Is vO.l solar masses. 

^, 7 ^- 119-19 
C.3/90 36995 
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I. Int roduot Inn 

Hni'auiii* IRl’ ♦102Jh In « prion* nxHmplt* of n o«rbi*n Htar undornolnj; 
oxtentilvc okiMM Iohh, a v.ood doal of ohaorvat lonal ami tlioorotloal atiidv 
ban ho#n Jovotod toward undoratandlntt tiu* atruoturo and dvn.imli'a of (ta 
c 1 ronmatollar onvolopo iMorrla, Kwan and Hill, l‘»77). Tin* Impor- 

lanoo of Bk'loonlar apool roaoopv In auoh a atndv la obvlona In that tin* 
a* 'noa of apootral llnoa art* dlroot Indlcatora of dvnamlcal aotlvlty. 

I'p to now tiu* boat voloolty roaolntlon liaa boon at raulo t ronnono loa , and 
oonaldorablo ovld. noo i>n nuaa li>aa liaa conu* frino obaorvatlona In thia 
apootral rottlon. Howovor, t ho ani;;ular ro'iolut Ion of aln>;lo radio t«*lo- 
aoopoa la onrront Iv Inaiif f lo lout ro roaolvo tho atruotnro within tho '1 
aromln nuiloonlar onvolopo; ami oonaoijuont Iv , tho 11m* pioflloa aro ob- 
aorvod to bo "v 10 km/aoc wldo, alnco thoy onoompaaa molociiloa ovor tho lull 
rango of projoctod oxpanalon volooltloa. A nwiro aonaltlvo probe of nu'loo- 
nlar donalty and radial voloolty varlatlona In tho oxp4imllii>; onvolopo 
would bo oxpootod from obaorvat Iona of Infrarod abaorptlon llnoa a>;alnat 
tho relatively anuill \1 aroaoo oontral oont Inuum aouroo. Tho abaorptlon 
lino profiles aro Influonood only bv nx'loouloa alonj; tiu* I lm*-ot -a I ^ht 
toward the central aouroo, ami thua thoy moatly almpllfv the Intorprot.i- 
t i on of the volooltv atruoturo. 

llio Introduction of hotorodvno apootroaoopy to tho Infrared now 
allows thoao abaorpt Ion llnoa to bo atudlod at a velocity roaolut Ion 
pi 'vloualv obtainable only at radio froquoncloa. I'hla l.ot tor reports tho 
Initial application of this technique to atollar apoct roacopy with tho 

detection of ammonia In tho c Ircumat ol lar onvolopo ot IRi’ +I02lb. Throe 

14 -I 

tranaltlona of NH^ In tho v^ rotat lon-vlbrat Ion band armind dSO cm 

have boon observed. Althou>;h amtiK'nla haa a well atudlod microwave lnvt*r- 

alon apoctrum near 24 iIHz and la froquontly found In Intoratollar clouds. 
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it ban so far escaped detection in stellar sources. The significance of 
its detection in IRC >10216 lies not only in adding one more molecule 
to the list, but especially in that aimnonia by itself can provide a 
unique probe for many parameters of Interest in the study of mass loss 
effects. 

II. Observations 

The observations were made on a number of nights in 1978 June with 
the McMath Solar Telescope of Kitt Peak National Observatory*. The ini- 
tial detections were obtained with a 81 cm auxiliary solar telescope and 
Improved upon a week later with the main 1.3 m telescope. At frequencies 
around 950 cm the beam widths are only slightly larger than the dif- 
fraction-limited sizes of 3.2 arcsec for the smaller telescope and 1.7 
arcsec for the larger. The heterodyne receiver employs a HgCdTe photo- 
diode mixer and a CO 2 laser local oscillator and is much the same as 
described by Betz et al. (1977). However, a new 64 channel filter bank 
of 20 MHz filters now extends the velocity coverage to a width of '^'14 km/scc. 
The system sensitivity with chopping and telescope losses is measured in 
each channel to be 1.0 x 10 for a signal-to-noise ratio of 1 after 
a 1 sec integration. Because the laser oscillates on only discrete rota- 
tion-vibration transitions of CO 2 and because only a limited intermediate 
frequency bandwidth of ''<1500 MHz is available from the mixer, only Doppler- 
sb'.fted transitions of NH^ in fairly close coincidence with available 
laser lines may be observed (Hillman et al. 1977). 


*Kitt Peak National Observatory is operated by the Association of 
Universities for Research in Astronomy, Inc., under contract with the 
National Science Foundation. 



In IRl* till* pitrt U'ulitr 11iu*m which were Ncarchcil tor wore I lio 
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rtKll,l>, .itul trrtUHllionh of Nil ^ In the v , 

rotrtt lon-vlhrat Ion hand around ^SO cm Kl)(ur<* I llluatratca I he oh> 
aerved prolllea for the tlrat three linea. The center velocity of each 
0,22 km/aec channel la accurate to *0.0S Wm aec. The a<Ht'«‘«^ line at 

cm ^ waa not detected to an equivalent width limit ot v 10 cm ^ 
In a bandwidth equal to that of the ohaerved linen. The hvpertlne com- 
pon^nta of ammonia are not reaolvahle In theae apectra and would not he 
dlacernlhle at hln*»*T reaolutlon, alnce the Peppier wldtha of the Indi- 
vidual componenta are verv much larger than the hvpertlne aplltttn>;a. 

To aaaure the correct Identification ot theae lines, each leatute 
waa ohaerved on several occasions over the course of IP Jaya. Phan^ea In 
the orhltal velocity of the Karth caused a trequonev shift ot '1 ch.innel 
per day for the line position, which contlrnu'd that the lines occurred In 
the correct sideband for ammonia. To Insure that no parts ot the wln>ta 
ot the NH^ lines are folded over from the opposite Infrared sideband 
after mlxln)t, all the line centers were observed at hl^h Intermediate 
frequencies '‘>PP MHr. It Is Important to emphasize that there are no 
baseline Irregularities In these spectra of the tvpe which commonlv 
plague microwave line receivers. The technical explanation cannot 
be qlven here. The continuum contribution from the opposite Infra- 
red sideband Is assumed to be uniform and has been c.iretullv m»‘asured 
and renu'ved In these sluttle-sldebaiul spectra. The conversion i^alns ot 
the two sldebandu are known to be equal to within •■O.St from laboratory 
measurements. Thus, apart from the noise level, the shapes ot the 
resolved, but Incouplcte, features are real. Arouitd the aijl^.-^ j«»d 
aO(b,h> line frequencies, wlilch are both In coincidence with ' 



UHIGINAL PAGE IS 
OK P(K)R QUALITY 


4 


laser lines, lunar spectra show only a flat continuum level with no 
spectral features from the terrestrial atmosphere. The aR(l,l) frequency, 
however, is 1440 MHz below the R(14) line of (the laser line). 

Absorption from the wide pressure-broadened wings of terrestrial CO 2 in 
2 alrmasses produces a amooth 50Z change in atmospheric transmission 
across the aR(l,l) spectrum. For this line only, a fully resolved 
calibration spectrum of the Moon was used to normalize the NH^ data 
from IRC -M0216. Because of Increased atmospheric attenuation and a shorter 
integration time, the noise level on this line is higher. Finally, there 
is some evidence that the rather narrow features marked A, B, and C may 
all be real. However, any confirmations will require further observations. 

III. Analysis 

A) H 2 Density and Rotational Temperature. 

The spectroscopy of ammonia is sufficiently complex that the 
population of rotational levels will probably not follow a simple thermal 
distribution in a circumstellar environment. Not only Is gas closer to 
the central star likely to be hotter, denser, and in a stronger radiation 
field than that further away, but the effectiveness of collisions in 
thermalizing the energy levels diminishes rapidly at larger radial dis- 
tances. Consequently, absorption lines from rotational levels requiring 
a high level of excitation may be formed at locations quite different 
from those where the low energy-level lines are produced. A thoroughly 
quantitative determination of temperature, density, and NH^ abundance as 
a function of radial distance can only be obtained from a careful analysis 
of resolved line profiles in combination with a refined model of mass dis- 
tribution in the circumstellar envelope. On the other hand, a reasonable 
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Ht.trt toward eHtlnuitlng thoae quant it Ich can bi* nuidt* from tho Iohh 
ambitious approach taken there. 

8 -3 

At low H., denMltlea, <10 cm , only the metUHtable rotational 

*■ 9 

levels (J“K) will be populated aa their radiative llfetlroea are "VlO aec 

(Oka et al. 1971). The non-metnatable levels (J>K) of ammonia can decay 

via allowed (AJ“-1, AK-0) tranaltlona in tlmea of 100 see or lea«. If 

radiative trapping is unimportant, the populations of metastable levels 

will then assume a non-Bolt zoumn distribution, and a unique rotational 

temperature will not be derivable. As the H., density is increased from 

this limit, more of the non-metastable levels can be populated via col- 

10 -3 

llslons; and at densities >10 cm the populations of levels up to 

about J~6 will all assume* a roughly thermal distribution. At inter- 

8 -3 10 —3 

mediate densities, 10 cm <p„ <10 cm , only some of the low energy non- 

a. 

metastable levels such as (2,1) and (3,2) will be coll Isionally populated, 
which will bring the population ratios between the metastable levels 
closer to the thermal limit. 

For IRC +10216, a reasonable upper limit for the hydrogen density 
near the aimaonia can be derived from the non-detection of the aQ(b,4) 
transition. The radiative lifetime of the (6,4) rotational level is 
about 1 sec, after which It decays via an 119 cm ^ photon to the (S,4) 
level. Since the aQ(6,6) transition was easily detected, and the exci- 
tation energy of the (6,4) level la only 25X higher, collisions must not 

be frequent enough to maintain any significant population In the (6,4) 

-15 2 

level. A collislonal cross section of MO cm for H., on NH^ and a 
kinetic temperature of 500 K (justified later) then implies that p. < 

”2 

lO^^cm”^, if radiative trapping at 119 cm ^ can be ignored. A lower 
bound on the H., density may be computed indirectly from infrared 
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DC4aurementa of the CO column density. If It may be assumed that the NH^ 

lies at or within the radius of the CO distribution. A CC/H 2 abundance 
-3 20 -2 

ratio of *^10 • a CO column density of 10 cm , and a radial dlstrlbu> 

15 8 *3 

tion over <10 cm then set a lower limit on the H., density of >10 cm 

A ^2 

(Barnes et al. 1977). Consequently, most of the NH^ is expected to be 

8 -3 10 -3 

observed at intermediate densities: 10 cm <p„ <10 cm . Actually, a 

”2 

slightly better estimate of the H 2 density can be obtained by comparing 

the relative Intensities of the 3 observed lines. 

Of the 3 line profiles of Figure 1, the most notable difference is 

between the aQ(2,2) and aQ(6,6) lineshapes. Although the aQ(6,6) is 

stronger than the aQ(2,2) in both the low and high frequency wings, it 

is weaker at line center. The aR(l,l) and aQ(2,2) profiles are more 

comparable, but still their relative intensities differ between the 

wings and the line centers. Since the aQ(6,6) and perhaps even the 

aQ(2,2) line appears optically thick at line center, a better estimate 

of the rotational temperature can be deduced from the relative strengths 

of the optically-thin wings between LSR velocities of -34 and -38 km/sec. 

Q ^2 10 —3 

For H 2 densities between 10 cm and 10 cm , the aQ(2,2)/aR(l,l) ratio 

by itself requires only that the rotational temperature, is >200 K. 

The aQ(6,6)/aQ(2,2) intensity ratio sets a better limit: 400 K < 

8 -3 

700 K. Furthermore, in the lloilt p„ <10 cm , the increased non-thermal 

"2 

population of the (2,2) metastable level would not allow the aQ(6,6)/aQ(2,2) 

ratio to be as high as that observed, even at much higher temperatures. 

8 -3 

The high ratio requires that p„ >10 cm and that the (3,2) and (2,1) 

“2 

rotational levels are also thermalized. Of course the above comparison 
makes the reasonable assumption that the abundances of the ortho- and 
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para- spec lea of aimonla are equal, and Ui.it tunh species can be charac- 
terized by the same rotational temperature. 

B) Location 

To date, all the moic.ules detected In the infrared have been asso- 
ciated with the Inner 2 arcsec of the clrcumstellar envelope (Hall and 
Ridgway, 1978). Most of the ammonia can also be localized In this region 
because of both the high rotational temperature and the high H., density. 
However, a more definite assessment of the gas distribution can be nuide 
after first considering the spatial distribution of the continuum source. 
The first measurements to resolve the Infrared continuum were the lunr.r 
occultatlon observations of Toombs et al.(1972). These results were inter- 
preted with a simple two component model: an opclcally-thick Inner shell 

of diameter 0.4 arcsec and temperature 600 K and an opt leal ly-th In shell 
(t 'vO. 2) of diameter 2 arcsec and temperature 375 K. Kach of these dis- 
crete shells was estimated to produce about half of the total continuum 
radiation at 1000 cm Recent Interferometric observations by Sutton 
(1978), however. Indicate that about lOX of the 900 cm ^ flux from 
IRC -fl0216 actually comes from an unresolved (<0.2 arcsec) component and 
that the remaining 90Z has a Gaussian Intensity distribution with a 1/e 
diameter of *v0.9 arcsec. These results Imply a more continuous process of 
mass loss rather than discrete outbursts. A similar conclusion has been 
reached by McCarthy (1978) and Selby and Wade (1978), who have measured 
the continuum distribution at a number of infrared wavelengths. Further- 
more, Crabtree and Martin (1979) have shown In a theoretical model that 

_2 

the R density distribution expected from a constant rate of mass loss Is 
also consistent with the original lunar occultatlon observations. 
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To entimate the ammonia dlatrlbutlon, a reasonable first approach 
might be to assume that the lines are ikeen In absorption against all of 
the 900 cm ^ continuum radiation. Unfortunately, this simple configuration 
leads to severe difficulties In Interpreting the optlcally-thlck appear- 
ance of the aQ(6,6) profile. The relatively weak saturation depth of this 
line would Imply that the vibrational temperature of the ammonia Is ’''96Z 
of the >400 K continuum brightness temperature. The upper level of the 
vibrational transition cannot be thermallzed at this temperature by colli- 
sions, as this would require H 2 densities In excess of the limit Imposed 
by the undetected aQ(6,4) line. The radiative lifetime of the upper level 
Is only 'V'O.l sec. Therefore, the population In the upper vibrational 
level would have to be maintained by the radiation field, which would 
approach the required Intensity only In the Immediate vicinity of the hot 
dust. However, If the ammonia were this close to the continuum source, 
then the observed line profiles would be many times broader than those 
observed. This is because outflowing gas over almost the entire 15 kn/sec 
range of projected expansion velocities would be seen In absorption 
against the finite size of the continuum source. Thus, most of the am- 
monia cannot be outside of the major part of the continuum source, and 
the Inclusion of some radiative trapping does not significantly alter 
this basic conclusion. This Interpretation Is strengthened by the obser- 
vation that the 'vX).9 fractional Intensity at the center of the aQ(6,6) 
line Is also the fraction of the total 900 cm ^ continuum radiation seen 
outside of 0.2 arcsec In the Interferometer results of Sutton (1978). 

90Z of the observed continuum Is therefore picked up from the 0.9 arcsec 
Gaussian continuum distribution while most of the detectable ammonia 
must be within this diameter. Now, If the vibrational temperature of the 
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amaonla wer# «v«n a« hot aa 600 K, any optlcally-thlck abaorptlon lint'a 
wo-jld be Intrlnalcally more than 80Z deep agalnat the even hotter Interior 
duat or the atellar photoaphere ar >1500 K. Thi Intenaity ratio between 
the apat tally reaolved and unreaulved continuum component a can alao be 
uaed to fix the optical depth of the duat in the reaolved part to be on 
the order of 1, if the unreaolved xource can be characterized at 'v>1500 K 
and the reaolved part at *v500 K. Optical deptha <1 are conaiatent with 
the obaerved Gauaalan diatrlbutlon for the reaolved continuum inter.alty. 

(At 119 cm Important for radiative excitation of the (6,4) level, the 
optical depth of the duat ahould be negligible.) 

Although moat of the aimaonia la Interior to a large part of the duat 
reaponaible for the 900 cm ^ continuum radiation, it ia atlll probably 
aaaociated with the region of duat formation. Thla la becauae radiation 
preaaure on duat giaina ia expected to drive the gas to the obaerved 
expanaion velocity, and no ammi'nla la aeen at H 2 denaitiea >10^^cm 
Tlie narrownesa of the profllea also requirea that the moleculca in the 
observed levela not be too close to the unreaolved (<0.2 arcaec) component. 
To some extent, this happens naturally, aa molecules closer to the star 
will be hotter and distributed over a larger number of energy levela. 
Preliminary calculationa also Indicate that, for such narrow lines, the 
ammonia abundance ahould at first increase with increasing diatance from 
the central source, as would be the case if the relative NH^ abundance 
were in equilibrium with the grain temperature and not fixed at the 
photospheric value. 

Finally, because the (2,2) and (1,1) rotational levels will have 
larger populations at low temperatures <200 K, the additional absorption 
seen in these lines below the 0.9 intensity level can be attributed to 
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opt Ically-thin NH^ f«r outMld« the warm dust producing the 400 cm 
continuum. In this respect It is Interest i||| to note the constancy of the 
expansion velocity in the asymmetric profile of the aQ(2,2) line. Most 
of the ammonia within the central 1 arcsec of IRC -4-10216 is observed at 
expansion velocities ^^tween 11 and 15 km/sec, the same range as other 
molecules found in the much larger "radio" envelope (Morris, 1975). Most 
of the acceleration in mass loss thus occurs within a few stellar radii. 
The remarkably low velocity dispersion seen in these lines (only a few 
times more than the thermal velocity of NH^ at 600 K) is in good agree- 
ment with the turbulence estimate derived by Geballe et al. (1973) based 
on observations of infrared uO lines. A good measurement of this turbu- 
lence is important for accurate modeling of the radiative transfer and 
cooling effects of various microwave and Infrared emission lines. 


C) Column Density 

After the uncertainties in rotational temperature, hydrogen density, 

and the possible saturation of the aQ(2,2) line are taken irto account, 

17 -2 

the radial column density of ammonia can be estimated at I x 10 cm 
from the equivalent widths of the aR(l,l) and aQ(2,2) lines, and the v., 
band intensity given by Taylor (1974). This estimate is uncertain by a 
factor of 2, and it assumes that only 10% of the observed continuum radia- 
tion is emitted behind the ammonia. 

17 -2 

The ammonia column density of 10 cm and the approximate volume 
8-3 15 

density of 'vlO cm at a radius of ‘vlO cm from the star can be used to 

24 -2 

get a rough estimate of the column density of 'vlO cm and a 

fractional abundance of -vio”^. This fraction is higher than the photo- 
spheric ratio calculated by Tsuji (1964) for a 1000 K supergiant, and 
may be further evidence that the ammonia abundance evolves after ejection 


frmi Wt\e pflbtoaph«ri.*. T1i« H. coluam dcnalty iapll<>ii tht* m.<HH uf 

• ® 

Chf clrcuiMCtfllar v^wlop* within a 1 arcaac radlua la about O.I aol/ir 

iLxaaa, an aatlnate which la Independent of any aaaumed abundance ratloa. 
It would be uaeful to coapare the relative abundancea of SH^ and CO. 
alnce the formation of both aoleculea la relatively Ineenaltlve to the 
amount of carbon loat to gralna. However, any real let Ic comparlaon 
would require a better underatandlng of the radial dlatrlbution of CO. 
a study which aav alao be poaalble with infrared heterodyne techniques. 

Thia reaearch waa supported in part by NASA Grants NCR 
and W7L 0)-003-272. We are grateful to A.C. Cheung and C.H. Townes for 
helpful diacuaalona, K.C. Sutton and U. Pltelson for technical assistance, 
and the staff of Kitt Peak National Observatory for their continued 


hospitality. 
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Figure 1 Caption 

14 

Sp«ctra of three NH^ lines in the \>2 band: aR(l,1) at 

971.8822 cm"^, aQ(2,2) at 931.3334 cn"^, and aQ(6,6) at 927.3232 cn"^ 
The Doppler shifts of these lines are plotted relative to the local 
standard of rest (LSR) . LSR velocities may be converted to heliocentric 
values by the relation: ■ V^SR ^ km/sec. The clrcumstellar 

expansion velocity is determined by te ing the '.ntridrlc stellar veloc- 
ity as -26 ± 1 kn/sec (Kulper et al. 1976). The observed intensities 
are expressed relative to the continuum level, and error bars are illus- 


trated for 2 standard deviations. 
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